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We calculate the transition temperature in 2 and 3-flavor QCD using improved gauge and staggered fermion 
actions on lattices with temporal extent N T = 4. We find T c — 173 (8) MeV and 154 (8) MeV for n/ = 2 and 3, 
respectively. In the case of 3-flavor QCD we present evidence that the chiral critical point, i.e. the second order 
endpoint of the line of first order chiral phase transitions, belongs to the universality class of the 3d Ising model. 



1. Introduction 

Determining the properties of strongly inter- 
acting matter at high temperature, the temper- 
ature at which the transition to the quark-gluon 
plasma phase takes place as well as other criti- 
cal parameters, is one of the basic goals in lat- 
tice studies of finite-T QCD. Numerical calcula- 
tions which have been performed during the last 
years have shown that the details of the transition 
strongly depend on the number of quark flavors 
(n/) as well as the value of e.g. the pseudo-scalar 
meson mass (mps) which is controlled by the bare 
quark masses (m q ). Furthermore, it became evi- 
dent that the moderate values of the lattice spac- 
ing (a ~ 0.25fm) used in finite-T calculations 
with dynamical fcrmions lead to sizeable cut-off 
effects; different discretization schemes used in 
the fermion sector, e.g. the standard staggered 
and Wilson fermion formulation, lead to signif- 
icantly different results for T c Calculations 
with improved actions thus seem to be manda- 
tory to obtain reliable quantitative results for T c 
and its dependence on n/ as well as mps- 

We will report here on results for nt = 2 and 3 
obtained with a Symanzik improved gauge action 
and a staggered fermion action with improved ro- 
tational invariance, i. e. the p4-action with fat 1- 
link terms. This action has been used previously 
by us in studies of the flavor dependence of the 
QCD thermodynamics Further details on the 
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action and the algorithm used for our simulations 
are given in || . For a discussion of recent results 
with improved Wilson fermions see Ref. ^ . 

2. Flavor and quark mass dependence of T c 

The transition temperature in QCD with dy- 
namical fermions has been found to be signif- 
icantly smaller than in the pure gauge sector. 
This is in accordance with intuitive pictures of the 
phase transition based for instance on the ther- 
modynamics of bag or percolation models. With 
decreasing m q hadrons become lighter and it be- 
comes easier to build up a sufficiently high par- 
ticle density that can trigger a phase transition. 
For the same reason such models also suggest that 
T c decreases when nj and in turn the number of 
light pseudo-scalar mesons increases. This gen- 
eral picture is confirmed by the lattice results pre- 
sented in the following. 

We have calculated the pseudo-critical cou- 
plings, f3 c (m q ), for the transition to the QCD 
plasma phase on lattices of size I6 3 x 4. Calcula- 
tions have been performed for 2 and 3-flavor QCD 
in a wide range of quark masses, 0.025 < m q < 
1.0. The smallest quark mass value corresponds 
to m PS ~ 350 MeV. 

In order to set a scale for the transition tem- 
perature we have performed calculations of the 
light meson spectrum and the string tension^] at 

2 We note that the heavy quark potential is no longer 
strictly confining in the presence of dynamical quarks. The 
definition of the string tension is based on potentials ex- 
tracted from Wilson loops which have been found to lead 
to a linear rising potential at least up to distances r ~ 2fm. 
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Figure 1. The critical temperature in units of y/a (left) and the vector meson mass (right) versus 
(mps/mv) 2 . Shown are results for 2, (2+1) and 3-flavor QCD obtained from calculations with the p4 
action on lattices with temporal extent N T = 4. For rif = 2 we also show results obtained by using 
unimproved gauge and staggered fermion actions. 



Pc(m q ) on lattices of size 16 4 . The results of this 
calculation are shown in Fig. 0. 

We note that T C J and T c /my do show a con- 
sistent flavor dependence of T c . In a wide range 
of quark mass values T c in 3-flavor QCD is about 
10% smaller than for rif — 2. The dependence 
on m q is, however, quite different when using the 
vector meson mass rather than the string tension 
to set the scale for T c . While T c /^fa does show 
the expected rise with increasing m q the contrary 
is the case for T c jmy. Of course, this does not 
come as a surprise. The vector meson mass used 
in Fig. ^ to set the scale is itself strongly quark 
mass dependent, my = m p + c m q . Its mass thus 
is significantly larger than the vector meson mass, 
m p , in the chiral limit. We stress this well known 
fact here because it makes evident that one has to 
be careful when discussing the dependence of T c 
on parameters of the QCD Lagrangian, e.g. rif 
and m q . One has to make sure that the observ- 
able used to set the scale for T c itself is not or 
at most only weakly dependent on the external 
parameters. 

In Fig. H we show results for the ratio \fa/mv 
calculated on 16 4 lattices for various values of the 
quark mass in 2, (2+1) and 3-flavor QCD. Also 
shown there is the result obtained in quenched 
QCD as well as from a partially quenched cal- 



culation which we have performed on gauge field 
configurations generated with a dynamical (sea) 
quark mass of m q = 0.1. The fact that this ratio 
shows only little flavor dependence and that the 
(partially) quenched limit is in reasonable agree- 
ment with experiment and phenomenology^ indi- 
cates that the corresponding partially quenched 
observables are suitable also for setting the scale 
in the presence of dynamical quarks. 

2.1. T c and thermodynamics of (2+1)- 
Flavor QCD 

Physically most relevant is an analysis of the 
thermodynamics of QCD with 2 light {m u ~ 
rrid ~ 0) and a heavier (strange) quark with m s ~ 
T c . Calculations of bulk thermodynamic observ- 
ables (pressure, energy density) suggest that the 
explicit quark mass dependence is small above T c . 
In particular, a recent calculation of the pressure 
with m s /T = 1 and m u ^d/T = 0.4 has shown 
that deviations from the case with three degener- 
ate quarks are only about 10% 0. This difference 
would be even smaller when m s rather than m s /T 
would have been kept fixed in the simulation. One 
thus can expect that already for T>2 T c bulk 
thermodynamics of the plasma phase will be well 
described by the equation of state of 3 degener- 

3 Systematic deviations are visible but are on the 10% level. 
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Figure 2. y/a/my versus (mps/my) 2 for 2, 2+1 
and 3-flavor QCD obtained with the p4-action. 
Also shown are results obtained with unimproved 
staggered fermions. The filled black squares show 
results from quenched and partially quenched cal- 
culations with a sea quark mass of m q = 0.1. 



ate, nearly massless flavors. This is different close 
to the transition which in the (2+l)-flavor case 
presumably is just a rapid crossover. Although 
this analysis has not yet been performed with a 
sufficiently light light quark sector the result for 
the pseudo-critical temperature shown in Fig. 1 
suggests that T c is close to that of 2-flavor QCD. 

2.2. Cut-off effects 

In Fig. 1 we also show results for T c obtained 
with the standard, unimproved gauge and stag- 
gered fermion actions. We note that on lattices 
with temporal extent N T — 4 this action leads to 
systematically smaller values for T c j %fa as well 
as T c /mv- This provides an estimate for sys- 
tematic cut-off effects, which for a wide range of 
quark mass values seem to be about 10%. Such 
an estimate also is consistent with the observa- 
tion that for large values of N T calculations with 
unimproved fermion actions do lead to larger crit- 
ical temperatures. Statistical errors are, however, 
large in this case. We thus expect that our cur- 
rent analysis with improved gauge and staggered 
fermion actions is affected with systematic cut- 
off effects which are significantly below the 10% 



level. However, this clearly has to be checked with 
simulations on larger lattices. 

3. T c in the chiral limit 

For rif — 3 the chiral phase transition is known 
to be first order |Q| whereas it is most likely a 
continuous transition for rif = 2. This also im- 
plies that the dependence of the (pseudo)-critical 
temperature on the quark mass will differ in both 
cases. Asymptotically, i.e. to leading order in the 
quark mass one expects to find, 

T c (m PS ) - T c (0) ~ m PS ' n f ~ 2 , (1) 

I "Jps , 71/ > 3 

with 1/(38 = 0.55 if the 2-flavor transition in- 
deed belongs to the universality class of 3d, 0(4) 
symmetric spin models. Our estimates of T c in 
the chiral limit are based on the data shown 
in Fig. 1. For rif — 3 we have extrapolated 
Tc/a/ct and T c /mv using an ansatz quadratic in 
mps/my. In addition we have extrapolated my 
calculated at (3 c (m q ) for m q = 0.025 and 0.05 
linearly in m q to the critical point in the chiral 
limit, P e (0) = 3.258(4). The extrapolation to the 
chiral limit is less straightforward in the case of 
rif = 2. The data shown in Fig. 1 indicate that 
the quark mass dependence for rif = 2 and 3 is 
quite similar. This suggests that sub-leading cor- 
rections, quadratic in mps/mv, should be taken 
into account in addition to the leading behavior 
expected from universality. In our extrapolations 
for n/ = 2 we thus also add a term quadratic in 
mps/mv to the leading term given in Eq. 1. 
From these fits we estimate 
Tc_ _ / 0.225 ± 0.010 , n f = 2 
m p ~\0.20±0.01 , n/ = 3 ' (2) 

which corresponds to T c = (173 ± 8) MeV and 
(154 ±8) MeV for rif = 2 and 3, respectively. For 
rif = 2 this agrees well with results obtained in 
a calculation with improved Wilson fermions || . 
We stress, however, that the errors given here as 
well as in Ref. 3 are statistical only. Systematic 
errors resulting from remaining cut-off effects and 
from the ansatz used for extrapolating to the chi- 
ral limit are expected to be of similar magnitude. 
In order to control these errors calculations on lat- 
tices with larger temporal extent are still needed. 
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Figure 3. Contour plot for the joint probability distribution of Sg and iptp (left) as well as E and M 
(middle) for 3-flavor QCD and the 3d, 3-state Potts model (right). The QCD contour plots are based 
on calculations performed on a 16 3 x 4 lattice at j3 = 5.1499 with m q = 0.035. These parameters are in 
agreement with a previous estimate for the location of the chiral critical point || . The mixing parameters 
have been fixed to r = 0.55 and s — 0. Parameters for the Potts model are given in Ref. 7. 



4. Universality at the chiral critical point 
of 3-flavor QCD 

The first order chiral transition of 3-flavor QCD 
will persist to be first order for m q > up to 
a critical value of the quark mass. At this chi- 
ral critical point the transition will be second or- 
der and it has been conjectured that it belongs 
to the universality class of the 3d Ising model 
Q. From a simulation with standard, i.e. unim- 
proved, gauge and staggered fermion actions we 
find support for this conjecture. 

The analysis of the critical behavior at the 2 nd 
order endpoint of a line of 1 st order phase transi- 
tions requires the correct identification of energy- 
like and ordering-field (magnetic) directions |||7| . 
The ordering field at the critical point can be con- 
structed from a linear combination of the gluonic 
action Sg and the chiral condensate tjjip, 



E = Sr 



M = 4)il> + s S G 



(3) 



Here the mixing parameter r is determined from 
the m g -dependence of the line of first order tran- 
sitions, r = (d/3 1 'dm q )~* dpoint and s by demand- 
ing (E ■ M) — 0. In fact, unlike for energy-like 
observables, e.g. critical amplitudes that involve 
the thermal exponent yt, the universal properties 
of observables related to M do not depend on the 
correct choice of the mixing parameter s as long 
as the magnetic exponent y^ is larger than y t . 
The joint probability distribution for E and M 
characterizes the symmetry at the critical end- 



point and its universality class. Contour plots for 
the normalized E-M distributions at the critical 
endpoints of 3-flavor QCD and the 3d, 3-state 
Potts model as well as the corresponding plot for 
the SG'iiJtp) distribution are given in Fig. 3. This 
shows that also in the QCD case a proper defini- 
tion of energy-like and ordering-field operators is 
needed to reveal the symmetry properties of the 
chiral critical point. The joint distributions of E 
and M suggest the universal structure of the E- 
M probability distribution of the 3-d Ising model 
|j| , although it is apparent that an analysis of 3- 
flavor QCD on larger lattices is needed to clearly 
see the "two wings" characteristic for the 3d Ising 
distribution. 
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